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Abstract 

Background: The associations between vitamin D receptor (VDR) gene polymorphisms and breast cancer risk were 
comprehensively investigated to clarify issues that remain controversial. 

Methodology/Principal Findings: An electronic search was conducted of several databases, including PubMed, the 
Cochrane library, Web of Science, EMBASE, CBM and CNKI, for papers that describe the association between Fokl, poly-A 
repeat, Bsml, Taql or Apal polymorphisms of the VDR gene and breast cancer risk. Summary odds ratios and 95% 
confidence intervals (CI) were estimated based on a fixed-effect model (FEM) or random-effect model (REM), depending on 
the absence or presence of significant heterogeneity. A total of 39 studies met the inclusion criteria. A meta-analysis of high- 
quality studies showed that the Fokl polymorphism of the VDR gene was associated with an increased risk of breast cancer 
(ff vs. Ff+FF, OR: 1.09, 95%CI: 1.02 to 1.16, p = 0.007). No significant associations were observed between the other 
polymorphisms and breast cancer risk. No positive results were detected by pooling the results of all relevant studies. 

Conclusion: A meta-analysis of high-quality studies demonstrated that the Fokl polymorphism of the VDR gene was closely 
associated with breast cancer risk. 
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Introduction 

Laboratories investigations and epidemiological studies have 
suggested that the level of vitamin D, and the expression of the 
vitamin D receptor (VDR), might be associated with an increased 
risk of breast cancer [1,2]. However, based on the current 
available data, these relationships need to be further evaluated. 
Vitamin D from all sources undergoes hydroxylation in the liver to 
become 25-hydroxyvitamin D [25(OH) D], which is then further 
hydroxylated in the kidneys and other tissues to an active form of 
vitamin D (1,25-dihydroxyvitamin D, l,25(OH) 2 D) [3,4]. In 
several studies, l,25(OH) 2 D has been demonstrated to promote 
cell differentiation and inhibit cell proliferation, potentially 
modifying cancer risk via binding to the VDR [5,6]. The VDR 
is an intracellular hormone receptor that specifically binds to 
l,25(OH) 2 D and interacts with specific nucleotide sequences 
(response elements) of target genes to produce a variety of 
biological effects. As vitamin D exerts its activity by binding to the 
VDR, the finding that normal breast epithelial cells [7] and most 
breast cancer cells [8] express VDR suggests the possibility that 
VDR gene polymorphism may be associated with breast cancer 
risk. 



The gene that encodes VDR maps to the long arm of 
chromosome 12 (12ql2-14), and harbors approximately 200 
single nucleotide polymorphisms (SNPs). Some are linked to 
differences in l-25(OH) 2 D uptake and can therefore be considered 
as latent disease risk variants. A series of characterized VDR gene 
polymorphisms, including Fokl (rs2228570) [9-24], a poly- 
adenosine (poly-A) repeat variant [10,19,22,25-27], Bsml 
(rsl544410) [10,11,13,15-17,19-22,24,25,28-33], Taql 
(rs731236) [9,12-14,17,18,24,26,29,30,34-38] and Apal 
(rs7975232) [9,13,17,18,24,29,36,38-40], have been extensively 
studied with regard to their association with breast cancer risk, but 
with conflicting results. To clarify the association between breast 
cancer risk and VDR gene polymorphisms, we performed a meta- 
analysis of 39 existing studies to clarify the relationship between 
genetic variations in VDR and the risk of breast cancer. 

Materials and Methods 

Search Strategies 

A comprehensive literature search of numerous databases, 
including PubMed, the Cochrane library, Web of Science, 
EMBASE, CBM (China Biology Medicine) and CNKI (China 
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National Knowledge Infrastructure), was conducted up until 
December 21 st , 2013. Publications with the following search 
words in the titles, abstract or key words of the original studies 
were included: 'vitamin D receptor', 'VDR', 'Fokl', 'Poly A', 
'Bsml', 'Taql', 'ApaV, 'polymorphism' or 'variant' or 'mutation' 
coupled with the term 'breast cancer'. Additional studies that were 
not captured by the database search were identified by reviewing 
the bibliographies of relevant articles. 

Inclusion Criteria 

All identified studies were reviewed independently by two 
investigators. The following criteria were used for a publication to 
be included in the meta-analysis: (1) any study published as an 
original study that evaluated the association between VDR gene 
polymorphisms [Fokl, poly A, BsmX, Taql and ApaV] and breast 
cancer risk; (2) cases of breast cancer were confirmed by medical 
records or linkage with population-based tumor registries; (3) the 
numbers of case and control groups for each genotype were 
reported or the relevant data was available, and adequate data was 
provided to calculate the odds ratio (OR); and (4) publications in 
both English and Chinese were included. 

Data Extraction and Quality Assessment 

Two investigators conducted the search, extracted and tabulat- 
ed all the relevant data independendy. If a study was referenced 
more than once, the most complete and newly released study was 
used. If one article reported two or more different case-control 
studies, it was considered as two or more studies, respectively. 
Data extracted from each study were as followings: name of the 
first author, publication year, ethnic origin of the studied 
population, numbers of case and controls, and the genotype 
frequency of the polymorphisms. To maintain consistency with the 
previously published literature, SNPs of the VDR gene were 
reported using restriction fragment length polymorphism (RFLP) 
nomenclature for the major and minor alleles, as follows: Fokl 
(rs2228570) alleles C = F and T = f; Bsml (rsl544410) alleles G = b 
and A = B; Taql (rs731236) alleles T = T, and C = t; and Apal 
(rs7975232) A = A and C = a. The allele counts were calculated 
from the genotype counts when needed. 

The quality of studies was assessed according to the STrength- 
ening the REporting of Genetic Association Studies (STREGA) 
criteria [41], and studies according with STREGA criteria were 
defined as high-quality studies. An independent review and 
decision was made by a senior investigator if there were 
disagreements between the two initial reviewers. 

Statistical Analysis 

The strengths of the associations between five polymorphisms of 
the VDR gene and the risk of breast cancer were assessed for the 
contrast between two groups of homozygotes ff vs. FF, SS vs. LL, 
bb vs. BB, « vs. TT, aa vs. AA), the recessive iff vs. Ff+FF, SS vs. SL+ 
LL, bb vs. Bb+BB, ft vs. Tt+TT, aa vs. Aa+AA), dominant ff+Ffvs. 
FF, SS+SL vs. LL, bb+Bb vs. BB, tt+Tt vs. TT, aa+Aa vs. AA) and 
allelic if vs. F, S vs L, b vs. B, t vs. T, a vs. A) models by calculating 
the pooled OR and its 95% confidence interval (CI). The pooled 
ORs were obtained using either the fixed-effects (Mantel- 
Haenszel' method) model [42] or the random-effect (DerSimonian 
and Laird method) model [43], depending on the absence or 
presence of significant heterogeneity. The significance of pooled 
ORs was determined by the Z test. Heterogeneity among studies 
was assessed by the Chi-square test -based Q statistic and was 
quantified using the f statistic [44] . A significant Q_ statistic (P- 
value <0.10) or I 2 statistic (I >50%) indicated significant 
heterogeneity existed across studies. 



A Sensitivity analysis was performed to evaluate the key studies 
that had substantial impacts on between-study heterogeneity levels 
by removing the individual studies sequentially. To further explore 
the cause of heterogeneity, a meta-regression was performed, 
which included covariates such as ethnicity and sample size of the 
studies. If the origin of heterogeneity was found, subgroup analyses 
were conducted according to the origin. All statistical analyses, 
except the meta-regression, were performed using RevMan 
version 5. 1.6 software (Review Manager, Copenhagen: the Nordic 
Cochrane Centre, The Cochrane Collaboration, 201 1). The meta- 
regression procedure was conducted using STATA statistical 
software (version 12.0; Stata Corporation, College Station, USA). 

The possibility of publication bias was assessed using Begger's 
linear regression and funnel plots. An asymmetrical funnel plot 
suggested a possible publication bias. 

Results 

Baseline Characteristics 

A flow chart of the literature search is shown in Figure 1. 
According to the criteria eligibility, 39 studies was identified 
regarding the associations between the Fokl, poly-A, Bsml, Taql or 
Apal polymorphisms of VDR gene and breast cancer risk. Among 
these studies, 22 studies [9-24] concerned the association of the 
Fokl polymorphism with breast cancer, including 16,353 cases and 
21,881 controls, while seven studies [10,19,22,25-27] investigated 
the association between the poly-A repeat variation and breast 
cancer risk, with 5,493 cases and 5,566 controls. For the Bsml 
polymorphism, 25 studies [10,11,13,15-17,19-22,24,25,28-33] 
included 16,160 cases and 21,023 controls, while 16 studies [9,12- 
14,17,18,24,26,29,30,34-38] on the Taql polymorphism included 
6,940 cases and 8,267 controls. For the Apal polymorphism, 1 1 
studies were included [9,13,17,18,26,29,36,39] with 3,738 cases 
and 4,489 controls. All of these 39 studies provided sufficient data 
to calculate the possible relationship between the five polymor- 
phisms of the VDR gene and breast cancer risk. The general 
characteristics of the selected studies are summarized in Table 1 . 
More detailed information is shown in Tables SI, S2, S3, S4, S5. 
The pooled results are shown in Table 2. 

Meta-analysis 

Fokl polymorphism and breast cancer risk. Sixteen of 22 
studies were in accordance with STREGA criteria and were 
therefore defined as high-quality studies [1 1-16,18,19,21,22]. The 
meta-analysis of these studies showed a significant effect of the ff 
genotype on risk of breast cancer ffvs. Ff+FF OR: 1.09, 95%CI: 
1.02 to 1.16, p = 0.007; I 2 = 18%, p h = 0.24) (Figure 2). No 
significant associations were found for the other comparisons (ff vs. 
FF OR: 1.10, 95% CI: 1.00 to 1.20, p = 0.06; ff+Ffvs. FF OR: 
1.03, 95% CI: 0.97 to 1.10, p = 0.33;/vs. F OR: 1.04, 95% CI: 
0.99 to 1.09, p = 0.12 (Figures SI, S2, S3). No positive results were 
detected by pooling the data from all 22 studies. 

Poly-A variant and breast cancer risk. Four of the seven 
studies complied with the STREGA criteria [22,25,27]. No 
significant association was detected between the poly-A variant 
and breast cancer by pooling the results of all studies or only high- 
quality studies. 

Bsml, Taql and Apal polymorphism and breast cancer 
risk. For all relevant studies, the pooled results did not illustrate 
any significant correlation between Bsml, Taql or Apal polymor- 
phisms and breast cancer risk. Twelve high-quality studies did not 
show a significant association between the Bsml polymorphism 
and breast cancer risk [11,13,15-17,19,21,24]. As less than three 
studies complied with the STREGA criteria for the Taql and Apal 
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relevant studies 
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35 duplicated records 
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4 review articles 
3 only include case group 
2 with insufficient data 
1 with repeated data 



Figure 1. Flow diagram of the literature search. 

doi:10.1371/journal.pone.0096125.g001 



polymorphisms, and therefore meta-analyses were not performed 
by pooling the results of high-quality studies alone. 

Meta-regression 

To detect the origin of study heterogeneity, the random effects 
meta-regression method was used [45]. In the regression 
procedure, an independent variable, the logarithm OR, and two 
covariates, ethnicity and sample size was included. The results of 
all the meta-regressions showed that the two covariates were not 
the origin of the heterogeneity. 

Sensitivity Analysis and Publication Bias Evaluation 

A sensitivity analysis was performed by removing the individual 
studies sequentially to assess the effect of individual studies. The 
results detected did not differ from the initial analysis. Begger's 
linear regression showed that no publication bias existed in 
relationship to any variation (P>0.05). The funnel plot for the 
recessive model of Fokl polymorphism was symmetrical (Figure 3). 
The other funnel plots were not shown. 

Discussion 

Vitamin D regulates a variety of independent biological 
processes including bone metabolism, the innate immune 
response, cell proliferation and cell differentiation [46,47] . Several 
studies have suggested that adequate vitamin D levels may provide 
protection against chronic diseases, such as cancers, and could 



improve cancer prognosis [48]. The important roles that VDR 
polymorphisms play in the pathogenesis of breast cancer have 
been investigated across the world. Many studies have been 
carried out to investigate the relationship between VDR gene 
polymorphisms and the risk of breast cancer [9-39]. As a result of 
the limitations of sample sizes and the low statistical power of 
individual studies, research results have been conflicting and 
inconsistent. A previous meta-analysis involved only four SNPs 
(Fokl, Bsml, Taql and Apal) and contained relevant studies that 
were published before October 2008 [49]. Another review 
published in 2009 summarized the association between VDR 
polymorphisms and breast cancer risk, but no definitive quanti- 
tative results were obtained [50]. Our current meta-analysis 
included almost all studies that had investigated the associations 
between Fokl, poly-A, Bsml, Taql and Apal polymorphisms and 
breast cancer risk. This is the most comprehensive meta-analysis to 
data, to the best of our knowledge. 

This meta-analysis included data from 39 relevant studies. The 
meta-analysis of high-quality studies showed that individuals with 
homozygous ff genotype were responsive to the increased risk of 
breast cancer compared to patients with Ff or FF genotypes. The 
overall data from all genetic models did not demonstrated that 
there was a significant association between the poly-A repeat, 
Bsml, Taql and Apal polymorphisms and breast cancer risk. A 
sensitivity analysis was performed by removing the individual 
studies sequentially, and the overall genetic effects were consistent 
with those of the corresponding sensitivity analyses for the poly-A, 
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Bsml, Taql and Apal variants. These findings further indicated the 
robustness of the lack of association between these four polymor- 
phisms and breast cancer risk. 

No linkage disequilibrium was shown between the Fokl 
polymorphism and Bsml, Taql and Apal polymorphisms [51]. 
Therefore, Fokl can be considered as an independent marker 
within the VDR gene. The effect of this SNP on breast cancer is 
plausible, given that the presence of the f allele in the 5'-promoter 
region of the VDR gene results in a protein that is three amino 
acids longer protein that the wild-type, and which is less 
transcriptionally active [52]. The presence of this polymorphism 
could therefore cause reduced effects of vitamin D. Although most 
previous studies on the association of Fokl polymorphism with 
breast cancer did not identify any evidence for a significant 
association, McKay et al. found a positive statistically significant 
association between ff genotype and the increased risk of breast 
cancer by pooling six studies (OR: 1.16, 95%GI: 1.04 to 1.28, 
p = 0.006) [16]. Similarly, a previous meta-analysis published in 
2009 also showed a significant increased risk of breast cancer in ff 
genotype carriers (OR: 1.15, 95% CI: 1.03 to 1.26, p = 0.010) [49]. 
A meta-analysis of high-quality studies, provided strong evidence 
that the ff genotype was significantly associated with risk of breast 
cancer, in accordance with previous reports [49]. 

The poly-A repeat in the 3 '-untranslated region of the VDR 
gene which is strongly linked with Bsml, Apal and Taql has an 
important impact on IT)R mRNA stability [52]. Several studies 
have indicated that LL genotype (long/long) confers susceptibility 
to breast cancer risk compared with the SS genotype [10,22,26]. 
Chakraborty et al. revealed that the LL genotype is significandy 
associated with high-grade breast cancer in northern Indians 
[(unadjusted OR (95% CI): 4.45(1.87, 10.63); adjusted OR (95% 
CI): 4.66 (1.88, 11.53)] [26]. However, this result conflicted with 
the report from Ingles et al., where breast cancer risk was found to 
increase with increasing numbers of S alleles [10]. Our finding did 
not show any significant association between the poly-A variation 
and breast cancer risk in any genetic model. These inconsistent 
results might result from differences of ethnicity, sample size, study 
design, amongst other factors. 

The Bsml polymorphism is located at the 3' end of the IDR 
gene. It does not appear to change the nature of the translated 
VDR protein [53]. However, this polymorphism is linked in a 
haplotype with the variable-length poly A sequence within the 3'- 
untranslated region, which affects the VDR mRNA stability [54]. 
On the other hand, the Bsml, Taql and Apal polymorphisms are 
all in the same linkage disequilibrium block. These polymorphisms 
have been widely investigated, but with differing results. Consis- 
tent with a previous meta-analysis, our finding showed no 
significant association of these three genetic variations with breast 
cancer risk. Several studies have been performed to examine the 
VDR haplotypes [9,13,14,18,35], but these results were also 
conflicting. McCullough et al. has investigated haplotypes that 
involved Bsml(B/b), Apal (A/ 'a), TaqKTIt) and a poly-A repeat(.S7 
L), but they failed to find significant association between any 
haplotype and breast cancer risk [13]. However, in a Caucasian 
population the baTL has been reported to increase the risk of 
breast cancer [9,35]. It is unclear whether chance or underlying 
differences in populations led to these inconsistencies. Due to the 
limited information available about these polymorphisms, we 
could not conduct an analysis for linkage disequilibrium and 
haplotypes. 

A few studies have investigated the association of VDR 
polymorphisms with breast cancer survival, but their results were 
also inconsistence. An analysis conducted among 111 Swedish 
breast cancer patients younger than 37 years of age found a trend 
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Figure 2. Forest plots of association of Fokl polymorphism with breast cancer. Significant association was detected between the genotype 
ff and breast cancer in recessive model (ffvs. Ff+FF). The squares and horizontal lines correspond to OR and 95% CI of specific study, and the area of 
squares reflects study weight. The diamond represents the pooled OR and 95% CI. Heterogeneity was checked by the chi square based Q test. 
doi:10.1371/journal.pone.0096125.g002 



towards a higher survival rate, especially among those estrogen 
receptor-positive tamoxifen-treated patients that were homozy- 
gous for the rare Taql allele [35] . However, Perna et al. reported 
that homozygous carriers of the rare Taql homozygous genotype 
had a 2.8-fold increase in the probability of death from breast 
cancer compared to homozygous carriers with the common allele 
(OR: 2.8, 95%CI: 1.1-7.2) [55]. 

Our meta-analysis illustrates strong evidence for the association 
between a VDR gene polymorphism in Fokl and an increased risk 
of breast cancer. The obvious evidence of between-study 
heterogeneity in this meta-analysis should be discussed. Although 
a meta-regression procedure that included two covariates was 
performed, the origin of the heterogeneity among the studies was 



not found. The heterogeneity might have been due to other 
factors, such as diversity in the population characteristics 
(ethnicity, age, sun exposure and dietary vitamin D intake, etc.), 
genotyping methods and study design. Previous studies have 
shown that the ethnic (genetic) background, gene-gene or gene- 
environment interactions and life-style (sun exposure, dietary 
vitamin D intake and smoking) might play a major role in the 
increased risk of breast cancer in association with genetic 
variations. Our meta-analysis was based on estimates without 
adjusting the data for these factors, which is another potential 
limitation of this study. 

In conclusion, this comprehensive meta-analysis of high-quality 
studies provides substantial evidence that the Fokl polymorphism 



Figure 3. Begg's funnel plot to examine publication bias for comparisons of Fokl polymorphism (ffvs. Ff+FF). 

doi:10.1371/journal.pone.0096125.g003 
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in the VDR gene is significantly associated with an increased risk of 
developing breast cancer. Furthermore, individuals that were 
homozygous for the minor allele genotype of Fokl were more likely 
to develop breast cancer. No correlations were found between the 
poly-A variation, Bsml, Taql and Apal polymorphisms in the VDR 
gene and the risk of breast cancer in this study. 

Supporting Information 

Figure SI Forest plots of association of Fokl polymor- 
phism with breast cancer (ff vs. FF). 

(TIF) 

Figure S2 Forest plots of association of Fokl polymor- 
phism with breast cancer (ff+Ff vs. FF). 

(TIF) 

Figure S3 Forest plots of association of Fokl polymor- 
phism with breast cancer (f vs. F). 
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